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ABSTRACT 


aneueenents of forward scattered light mene ~used to 
determine the particulate behavior in the combustor and 
across the exhaust nozzle of a small solid propellant rocket 
MOeOr . A MALVERN 2600 was successful in measuring 
multimodal particle distributions, and a locally designed 
system was successfully applied to monomodal distributions. 
The propellant tested was a reduced smoke propellant 
containing one-percent zirconium carbide by weight. The 
results suggested that ZrC reacts to completion on or very 
near the propellant surface. Zirconium oxide (Zr0>) 
agglomerates are believed to form on the surface. Surface 
agglomerate size decreased with increasing pressures to 
approximately 500 psig. Agglomerates broke up in the nozzle 
convergence to approximately 18 microns, independent of the 


pressure or nozzle entrance particle size.- Particle 
breakup occurred in the exhaust plume. System limitations, 


sources of error and suggested improvements are discussed. 
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enc ODUC ELON 


The use of metallized propellants offers the advantage 
of improved solid rocket motor performance in terms of 
increased energy, specific impulse, and suppression of high 
frequency acoustic combustion instabilities. Aluminum is 
currently the primary metallic additive of choice when 
exhaust plume visibility is not a mission consideration. 
However, the requirement for reduced smoke propellants for 
use in tactical air-to-air applications has dictated the use 
of alternate high energy propellants which utilize only 
small amounts of additives such as zirconium, zirconium 
carbide, and aluminum oxide to provide acoustic 
stabilization. 

The improvements achieved via such additives are not 
without penalty, certain loss mechanisms are introduced. 
The specific impulse efficiency of a propellant containing 
metal is generally less than that of a base propellant due 
to incomplete combustion of the metal in the motor and the 
presence of condensed metal and metal oxides in the nozzle 
flow. Particle-gas velocity lag represents the primary 
source of performance loss within the exhaust nozzle. 
Thermal lag generally has less effect on performance. Yet, 
neglecting thermal and velocity lag as loss mechanisms, the 


simple fact that the condensed particulates do not expand as 


the other combustion gases expand, results in a reduction of 
available specific impulse. 

Netzer et al. [Ref. 1], point out that the behavior of 
metal as it leaves the burning surface is not a well-defined 
phenomenon. Some particles upon reaching the burning 
surface depart immediately, while others agglomerate on the 
surface before passing into the gas flow. The mechanics of 
metal droplet combustion as discussed by Schorr [Ref. 2], 
indicate that resultant particles in the range of 1-2 
microns and smaller are of significance in determining 
exhaust plume signature, while those particulates of 5 
microns and greater are important in the determination of 
two-phase flow losses. The problem of predicting 
performance depends on an accurate determination of these 
two-phase flow losses which are a strong function of 
particle size. The aluminum oxide that results from the gas 
phase combustion of aluminum is generally less than two 
microns in size. These small particles have little velocity 
lag. Thus, making them smaller would serve little benefit. 

Several complex computer codes attempt to model the 
combustion process in order to predict solid rocket motor 
performance. The SPP (Solid Performance Program) ([(Refs. 
3,4], 1s one such code. These programs are semi-empirical 
in nature and require an a priori knowledge of particle size 


distribution entering the exhaust nozzle. This distribution 


must be experimentally determined, as the problem is 
currently intractable with regard to analytic methods. 

The following experimental methods have been employed in 
previous work at the Naval Postgraduate School: 

(a) collection of exhaust products 

(b) light transmission measurements 

(c) light scattering measurements 

(d) holography and high speed motion pictures. 
A summary of previous research conducted at the Naval 
Postgraduate School combustion laboratory is contained in 
Reference 1. The present investigation focused on the use 
of light scattering techniques to determine particle size. 
The method is non-intrusive and suited to measurements in 
both the motor and exhaust nozzle environments. The speed 
of the flow precludes successful application of most 
holographic techniques within the exhaust nozzle. 

Heretofore most studies have been directed at aluminized 
composites, and much is known about the combustion of 
aluminum. The behavior of other metallic additives have not 
been as fully explored. For example, zirconium carbide is 
used as a stability additive. However, it is not known 
whether surface agglomeration occurs, or whether the 
zirconium burns to zirconium oxide. Both of these processes 
will affect the size distribution and density of the 


particles in the gas phase. These in turn affect the 


ability of the particles to provide damping for a specific 
frequency of pressure oscillation. 
The goals of this investigation were: 


(1) to refine previously developed experimental and 
diagnostic techniques, and 


(2) to use the apparatus to determine the effects 
operating pressure and position within the motor and 
exhaust nozzle on the particle size distributions 
from a reduced smoke propellant. 

The propellant studied was a low smoke, composite which 
contained one percent Zirconium carbide by weight. The 
behavior of particulates at various positions in the motor 


and exhaust nozzle was examined as a function of chamber 


pressure. 


IIT. THEORETICAL BACKGROUND 


The performance of a solid rocket motor utilizing 
metallized propellants is significantly dependent on the 
diameter of the condensed particulates present in the flow. 
A non-intrusive, in-situ method is the most desirable means 
to determine particle size, since the measurement is real 
time and the environmental conditions (pressure, 
temperature, gas composition, etc.) are undisturbed. The 
method chosen for the present research was the diffractively 
scattered light technique. Very large particles cannot be 
measured with this technique due to the transmitted light 
interference with the light which is scattered at very small 
forward angles. Very small particles scatter at larger 
angles, but the intensity is low, making detection 
difficult. Absorbtion and refraction also become important 
for very small particles. Particle concentrations are also 
limited to those which do not exhibit significant multiple 
scattering effects. Within these limits, the mean particle 
size is readily determined from diffraction techniques, but 
size distributions cannot easily be determined. In general, 
a distribution type must be selected a priori, i.e., Rosin- 
Rammler, log-normal or other such distribution must be 
assumed. Yet, the diffraction technique is non-intrusive 


and can be~- used in the internal motor and nozzle 


environment, where velocities and particle concentrations 
can be high. 

The scattering and absorbtion of a plane monochromatic 
electromagnetic wave by an isotropic homogeneous sphere is 
described by Mie's solution to Maxwell's equations (Ref. 5], 
and the general theory of scattering developed by Mie is 
presented by Van de Hulst [Ref. 6]. According to Mie's 
theory, the characteristics of scattered light depend upon 
the particle size, shape, internal structure, refractive 
index with respect to the surrounding medium, the scattering 
angle and the incident wave polarization and wavelength 
(Ref. 7]. The solution to Maxwell's equations for multiple 
scattering in a concentrated sample, where interaction of 
scattering centers is significant, presently represents an 
intractable problem, even if boundary conditions could be 
specified. Even the solution without multiple scattering is 
quite complex and difficult to incorporate into diagnostic 
equipment. 

When the particles are much larger or much smaller than 
the wavelength of the incident light, the equations simplify 
to Fraunhofer diffraction and Rayleigh scattering, 
respectively. The light scattered from particles within 
these restricted ranges can be measured and used to readily 
determine particle size. The size of particles of interest 
in solid propellant combustion depend on the application. 


When considering two phase flow losses, particle sizes of 


interest are generally five microns or greater. Jones [Ref. 
8], has shown that Fraunhofer diffraction can be used for 
particles as small as one micron when a He-Ne laser is used. 
Fraunhofer diffraction is applicable in the case of larger 
particles because the forward scattered light in ee center 
lobe is due principally to diffraction. The mechanisms of 
reflection and absorbtion do not significantly contribute to 
the center lobe in this instance (this was’ shown by 
Hodkinson [{Ref. 9]). If these latter effects are present to 
a Significant extent, then the complete Mie equations must 
be used. 

The classical results of Fraunhofer diffraction from 


particles, i.e., the discrete rings associated with a 


monodispersion, are typically not observed in practical 


combustion experiments. Typically, the particulates are 
polydispersed amg, in many cases, a multimodal 
polydispersion is present. The intensity distribution due 


to a monodispersion is given by Fraunhofer diffraction [Ref. 
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where J, is the first order Bessel function, 4 =7™D/), 0 is 
the scattering angle for the forward direction, D is the 


particle diameter, and } is the wavelength of the incident 


Lake pee - When considering small angles (8 < 15%), the 


expression Simplifies to, 
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Dobbins [Ref. 11] showed that the normalized integrated 
intensity of forward scattered light from a polydispersion 


of large particles is given as, 
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where N;(D) is a distribution function describing the 
proportion of particles with diameter D in the _ sample. 
Dobbins et al. [{Ref. 11] used the Upper Limit Distribution 
Function (ULDF) developed by Mugele and Evans [(Ref. 12] and 
the same approach was followed in a portion of this 
investigation. The volume-to-surface mean diameter (Sauter 


mean diameter) D35 of a polydispersion is defined by, 


3 
ii N(D)D aD 
"32 15 a (4) 
af N(D) D°aD 


Dhe Ssigmiticance of D5 “versus Ds3, etc.) is that the 
illumination profiles for different size distributions of 
the ULDF plotted against (17D396/\) are brought into near 
meine 1 detice (Ref. 11}. 

Buchele [Ref. 10] showed that for 7D39¢/i less than 
three, a gaussian curve can be used for the scattered light 
intensity which closely matches the theoretical intensity 
profile for a polydispersion obtained by integrating 


Equation (3). The gaussian expression presented by Buchele 


is given as, 





= exp[-(.571D359/,) 2] (5) 


The data reduction technique developed by Netzer and Harris 
(Refs. 1,13] is based on Buchele's' work. The method 
involves comparison of the experimentally obtained intensity 
profile with that produced by Equation (5). Use of Equation 
(5) requires the value of I1(0), the forward scattered 
intensity at theta equal zero. Since transmitted light 
dominates at very low scattering angles, I1(0) cannot be 
measured. A ratioing technique eliminates this problem. 
Equation (5) can be used to ratio I(6))/I(6,) as follows, 


cee) 2 57) 2 Z 
2 = exp(-D32 (4-9 2( 63 - 8])) . (6) 








This eliminates the need to determine I(0). However, I( 61) 
at the "minimum" forward angle ©, must be taken from the 
experimentally measured profile at a location where the 
Peanenieten light effects are not significantly present. 
The resultant profile I(95)/I(°1) versus % can then be 
generated for a specified D359. D359 1S varied interactively 
until the profile generated by Equation (6) matches the 
measured profile. This procedure yields D35. However, this 
approach has an obvious disadvantage in that it is only 
accurate if the polydispersion is monomodal. 

A Malvern 2600 particle sizer was also used in the 
present investigation. The Malvern software is similarly 
based upon Fraunhofer diffraction. However, it is based on 
the inversion of the measured scattered data to obtain a 
particle size distribution. The Malvern allows use of a 
selected type of particle size distribution (Rosin Rammler, 
log~normal, normal) or of a model independent mode. The 
model independent method is described by the Malvern 2600 
Reference Manual [Ref. 14] as a constrained non-linear least 
squares method, which is able to measure multi-modal 
particle distributions with high resolution. The obvious 
advantages of the Malvern 2600 are the ability to measure 
multi-modal particle size distributions, and it has been 
independently tested. Dodge et al. [Ref. 15] conducted a 
comparison of the Malvern device and a phase doppler system 


and found that there was excellent agreement in point 


10 


measured Sauter mean diameters between the deconvoluted 
laser diffraction data and the phase doppler data. The 
Malvern system is currently one of the most widely used 
laser Aer raction devices andmteschows gneat promise for 
application to particle sizing in solid rocket motors. 

Other work has been done in the area of analytic methods 
that appears to provide the basis for refinements in 
Baitsc le Sizing and determination of particle size 
distributions. Koo and Hirleman [Ref. 16] evaluated the use 
of integral transform techniques on Fraunhofer diffraction 


Paceerns. These techniques offer the advantage of direct 


inversion of scattering data to obtain a size distribution. 


ia 


III. EXPERIMENTAL APPARATUS 


A. INTRODUCTION 

Laser diffraction systems have become widely used in the 
field of particle and droplet size measurements. Fraunhofer 
diffraction is the underlying principle involved. Hodkinson 
(Ref. 9] demonstrated that for particles greater than 
several microns in size, the primary contribution in the 
near forward small angle region is from diffraction rather 
than reflection or rerraction and Jones [Ref. 8] has proved 
that in this same region, Fraunhofer diffraction may be 
applied to particles as small as one micron when a helium- 
neon laser is chosen as the source of illumination. Figure 
3.1 details a simple schematic for a laser light diffraction 
particle sizer. Such a system typically consists of a He-Ne 
laser source, beam expander (spatial filter and collimating 
lens), Fourier transform lens, and a detector. 

Typically, small particles diffract light at large 
angles and large particles diffract light at small angles. 
Figure 3.2 depicts the diffraction process. Note that the 


undiffracted light is focused to a single point. 


Bo ROCKET Meter 
The series of experiments conducted required the use of 
two small solid propellant rocket motors. The motor 


depicted in Figure 3.3 is the same as that used by Pruitt 
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Figure 3.3 Motor Components 





Figure 3.4 Short Motor Components 
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PRetee gale The chamber diameter was 2.00 inches and the 
propellant was cylindrical, 1.98 inches in diameter and 
1.00-1.25. inches thick. All tests conducted were done with 
end burning grains. The long motor (Fig. 3.3) was used to 
obtain in-situ measurements of particle size prior to the 
nozzle entrance and in the exhaust. Figure 3.4 depicts the 
short motor which was used to obtain in-situ particle size 
measurements just aft of the propellant burning surface. 
This particular motor had been used previously in connection 
with holographic imaging. For additional specifications see 
Reference 18. 

Both motors were windowed, and the windows were kept 
clean using a nitrogen purge. The nitrogen was "spread" 
using a sintered metal sleeve which was inserted into the 
aperture. 

Ignition was accomplished using a BKNO3 ignitor which 
was fired by means of a nichrome filament energized by a 12 
VDC source. 

The propellant was bonded into the motor casing with a 
self-vulcanizing silicone rubber (RTV). The propellant 
studied was a reduced smoke propellant with 87% solids 
loading. Table 3.1 presents the major propellant ingredient 


weight percentages. 


C. LASER LIGHT SCATTERING APPARATUS 
The locally designed laser diffraction system makes use 


of an 8 mWatt He-Ne source in the same general configuration 
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EES Ignis a 


PROPELLANT INGREDIENTS 


Ingredient Approx. Wt. % 
R45M Sie 
RDX 4.0 
AP 82.0 
Zee 1.0 
megrcted in Figure 3.1. (A complete description of the 


apparatus is contained in Reference 20 [see Figures. 3.5 and 
Bro ||) . The scattered light is passed through a wide-pass 
laser line filter and a 500 mm focal length Fourier trans- 
form lens, which focuses the light onto the detector. The 
detector is a linear photodiode array (Reticon RL 1024G) 
which has 1024 photodiodes on 25 micron centers. The array 
1s mounted on a "mother" board which contains the sample- 
and-hold circuits. Each photodiode measures the intensity 
of the scattered light and outputs a proportional voltage. 
The array 1s scanned at a rate of 33 kHz and a 6 msec blank- 
ing pulse is used to sequence the sweeps. The range of 
scattering angles within the field of view is determined by 
a number of factors. The upper limit of scattering angle is 
determined by the diameter and focal length of the focusing 
lens, the distance between the focusing lens and the 
particles, the diameter of the motor window, and the 


Hemene/posittion GL the diode array [Ref. 1). The range of 


iy: 


scattering angles for this investigation was .01-.06 
radians, corresponding to an offset of 5 mm of the array and 
.02-.07 radians corresponding to a 10 mm offset. Because of 
the very small detectors, the system is not capable of 
measuring very weak scattering intensities. 

The second system utilized was a commercially available 
particle size analyzer manufactured by MALVERN Instruments. 
The particular instrument used in this series of experiments 
was the MALVERN 2600 HSD. Figure 3.7 depicts the MALVERN 
2600 system, which utilizes a 2 mW He-Ne laser with a 
collimated beam diameter of 9 mm (Table 3.2 details laser 
specifications). One of the principle advantages of the 
MALVERN system is its use of an annular photodiode array. 
The array is a monolithic detector which consists of 31 
rings. It has increasing thicknesses of the annular detec- 
tor elements with increasing circumference, resulting in an 
increase of detector area as radius increases. Ring arrays 
permit much weaker intensities to be measured than can be 
accomplished with linear arrays. The dynamic range of the 
instrument is 180:1 on any of the three range settings, with 
an advertized accuracy of plus or minus 4% of volume mean 
diameter [{Ref. 19]. The detector may be used in conjunction 
with any one of three interchangeable Fourier transform 
lens, with focal lengths of 63, 100 and 300 mm [Ref. 14]. 
The respective size classes associated with each lens are 


1.2-118, 1.9-188 and 5.8-564 microns, with a sub-class 
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Figure 3.6 Hewlett Packard HP 9836S Computer and 


HP 6942A Multiprogrammer 
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LASER SPECIFICATIONS 


A. Helium-Neon Laser (Locally Designed System) 


1. Manufacturer: Spectra-Physics 
a Model: 147 

3. Type: He-ne Class 111B 
4. Output Power: 8 mWatt 

5. Beam Diameter: .92 mm 

6. Beam Divergence: .87 mrad 


B. Helium-Neon Laser (MALVERN System) 


1. Manufacturer: MALVERN Instruments 
2. Model: 2600 
ope: He-Ne Class 111B 
4. Output Power: 2 mWatt 
5. Beam Diameter: 9 mm 
6. Beam Divergence: Not Available 
which extends down to .5 microns [Ref. 14]. The 100 mm lens 


was used for this set of experiments. 


D. DATA ACQUISITION AND REDUCTION 

The data acquisition system associated with the locally 
designed laser diffraction system consisted of a Hewlett 
Packard 6942A multiprogrammer which was controlled by an HP 
9836S computer. The data acquisition program developed by 


Harris [{Ref. 13] was the system driver in the acquisition 
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mode. In the course of a motor firing sequence, a 
background measurement was recorded prior to ignition. This 
mean intensity profile was then subtracted from the mean 
intensity. profile measured with particles present. This 
procedure accounts for noise present due to ambient light 
and corrects for the characteristics of individual diodes. 
The raw data were displayed on the HP 9836S CRT so that any 
erroneous scans could be excluded prior to data reduction. 

A symmetric moving average type digital filter was then 
applied to the profile in order to smooth the curve [Ref. 
ais This type filter does not have the phase lag that is 
inherent in analog filters. The data reduction algorithm is 
based on the work done by Buchelle [{Ref. 10], and is fully 
described in Reference 20. The chief disadvantage of the 
method lies in the fact that there is no way to account for 
the actual distribution of the particulates. The Sauter 
mean diameter (SMD) determined by this algorithm will be 
that which 1S associated with the dominant mode (within the 
field of view of the array) of a multimodal polydispersed 
sample. 

The MALVERN system represents a vast improvement, in 
that not only is SMD determined, but the volume distribution 
is characterized as well. The MALVERN 2600 data acquisition 
system is centered about an AT&T 6300 PC. The system driver 
is the MALVERN Master-Sizer operating system. This version 


of the MALVERN software uses a command-based language that 
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allows the user to specify the pertinent parameters 
associated with the operation of the instrument. The 
Master-Sizer program allows the user to select from any one 
of four models. The model types available are Rosin- 
Rammler, log-normal, normal, and model independent. The 
model independent option of the Master-Sizer program is a 
process of constrained least squares fitting of theoretical 
scattering characteristics to observed data [Ref. 14}. Tete 
is usual to allow the best fit result to be obtained with no 
assumed form of size distribution and this allows the 
characterization of polydispersed multimodal distributions 
with high resolution [(Ref. 14]. The capability of 
characterizing multimodal distributions is invaluable in 
work dealing with the combustion of metallized propellants, 
in that the distribution of particulates 1s most often 
multimodal. It is for this reason that the model 
independent mode was used exclusively throughout this 


investigation. 
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Ly oe eeoUrls AND. DISCUSSION 


EX. INTRODUCTION 
This investigation was based on three major objectives: 

(i) “refine previously developed experimental and 
diagnostic techniques by identifying specific causal 
factors of measurement error, 

(2) incorporate the use of a commercial particle size 
analyzer (MALVERN 2600) to validate and quantify the 
operating limits of the locally designed laser 
diffraction system and, 

(3) use the apparatus to determine the effects of 
operating pressure and position within the motor and 
exhaust nozzle on Pause 1Ccle size Gist hole lon 
associated with the combustion of a reduced smoke 
propellant. 

Previous investigations at NPS had been based on the method 
of Buchele [{Ref. 10] which yields Sauter mean diameter, but 
does not provide any information as to the actual 
Gectribltaon of particulate size. mie sintweductaon of the 
MALVERN system allowed for a characterization of volume 
Custer 1 bution. The common parameter between the locally 
designed system and the MALVERN was Sauter mean diameter 
(D329). This parameter provided the basis for comparison of 
the two diffraction systems. Both systems were used in the 
initial evaluation of the reduced smoke propellant. Once 
the limitations of the locally designed instrument were 
established, the MALVERN was used more extensively to gather 


data in both the motor and exhaust. The distinct advantage 


of being able to characterize multimodal particle size 
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distributions made the MALVERN the better suited instrument 
for the through motor measurements. In the exhaust, the 
annular detector of the MALVERN system allowed for the 
measurement of much weaker intensities as opposed to the 
locally designed system (which is also referred to as the 


NPS system). 


B. SYSTEM CALIBRATION 

Among the methods available for calibration of laser 
@iffraction instruments, the more commonly used include the 
calibration reticle and use of commercially prepared samples 
of Known diameter. The latter approach was used in this 
investigation. The types of samples used included 
polystyrene spheres, glass beads and aluminum oxide powders. 
In each case, the sample material was suspended in a 
solution of distilled water. 

The locally designed system was calibrated using the 
polystyrene spheres, with diameters of 5.1, 9.6 and 20.0 
microns, respectively. Table 4.1 contains the results of 
these calibration tests. Figures 4.1 and 4.2. are 
representative samples of calibration results. 

The MALVERN 2600 was calibrated using three different 
materials, polystyrene spheres, glass beads and an aluminum 
oxide powder. Results are presented in Table 4.1. 
Calibrations from the two systems were in good agreement. 
System alignment was completed with the container of water 


in the optical path and a background measurement was’ taken 
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prior to introducing particles. A glass rod was used to 
stir the mixture to ensure that the particulates remained in 
suspension during the measurement process. 

A primary concern early in the investigation was to 
validate the ability oof the MALVERN to - accurately 
characterize a multimodal distribution when the system was 
operated in the Model Independent mode. In order to test 
this facet of the Master-Sizer program, a polydispersed 
bimodal sample consisting of 5.1 and 20.0 micron particles 
was prepared. The results are depicted as a volume 
Gistributionieim Figure s4-.. The respective modes are 
centered at approximately 5.5 and 21.0 microns. 

When multiple scattering becomes significant, the mean 
particle diameters determined from single scattering theory 
are too small. Gulder [Ref. 28] has determined empirical 
corrections for the Malvern when obscuration is greater than 


0.5. This correction was utilized, when required, in the 


present investigation. 


C. RESULTS FROM REDUCED SMOKE PROPELLANT EVALUATION 
1. Test Conditions 
The series of experiments conducted covered a range 
of pressures from approximately 70 to 450 psig. The 
experiments were divided into three pressure classes, 150, 


250 and 400 psig. At each class of chamber pressure, 
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scattering measurements were made at’ the following 
locations: 

(a) near grain (end burning) 

(b) betes to the nozzle entrance 

(c) exhaust (Pe = Peritical) 

(d) exhaust (P.e°= Po). 
A total of 32 tests were completed, 24 of which produced 
usable data. Three tests failed to produce data because the 
burst disk ruptured, three tests of the short motor were 
invalidated because the nitrogen purge system failed to keep 
the window surfaces clean, one test failed because of an O- 
ring failure in the motor joint, and for one test, data were 
not obtained because the pressure threshold to trigger the 
acquisition of data was set too high. 

The nozzles were designed/sized for isentropic flow 
conditions using y = 1.2. The ratio of specific heats was 
determined using PEPCODE, an equilibriun, adiabatic 
combustion code. The value of gamma cbdtained in this manner 
was 1.21. Because nozzle efficiencies should be quite high 
(97-99%), the isentropic assumption was made. Deviations in 
pressures can be attributed to the variations in ignition 
characteristics and propellant properties from run to run, 
and some small amount of nozzle throat erosion. 

The propellant composition was provided by the Naval 
Weapons Center at China Lake, California. The zirconium 


carbide in the propellant varied in size from 2-20 microns, 


Be 


with a D359 of 5.1 microns. The information is provided here 
as an INCrEOAUCE OI—eS the discussion of surface 
agglomeration on the burning grain. 

De Zirconium Carbide Behavior 

The size of agglomerates that form on the eamedee of 
burning metallized solid propellants are determined by many 
factors, including metallic powder type, propellant type, 
size of original metal particulates, concentration of metal, 
motor chamber pressure, and the burning rate of the 
propellant. The effect of chamber pressure on agglomerate 
size has been observed to be quite strong for aluminum. As 
chamber pressure and burning rate increase, the agglomerates 
become smaller [{Refs. 21,22]. 

Zirconium carbide has been observed to burn 
vigorously on the surface of reduced smoke propellants 
Containing AP and HTPB [Ref. 23]. If zirconium carbide 
(ZrC) burns stoichiometrically with oxygen it results ina 
temperature of approximately 4400°K. This resultant 
temperature is much greater than the melting temperature of 
ZrO> (2988°K). Table 4.2 delineates the temperatures 
relevant to the discussion. 

A plausible behavior of the Zrc is that as the 
material is exposed at the grain burning surface, oxygen 
(O>) from the ammonium perchlorate (AP) flame attacks the 
carbon to produce either carbon monoxide (CO) or carbon 


dioxide (C05). Oxygen from water vapor may also be 
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TABLE 4.2 


THERMOPHYSICAL PROPERTIES OF Zr, Z2rC, AND Zr0O5 


a Melting Boiling 
Material Temperature (°K) Temperature (°K) 
Zr 2127 [Ref. 24] 4700 [Ref. 24] 

Zxrc 3523 [Ref. 26] a 
Zr0> 2988 [Ref. 24] 4573 (Ref. 25] 


NOTE: Adiabatic flame temperature for the propellant was 
determined to be 2916°K and the adiabatic flame 
temperature for the reaction 

ZCC + 205 Be ZrO5 5 CO5 


is approximately 4400°K. 


involved. The oxygen then reacts with the zirconium to 
produce a porous ZrO, surface layer {Ref. 24]. The oxide 
should be a liquid since the local temperature at the 
particle surface (as a result of the surface combustion of 
ZrC or Zr with O5) should be much greater than the melting 
Eenperacure or ~ZEe> (2988°RK). The low burning rate 
associated with the relatively low chamber pressures allows 
time for the molten ZrO5 on adjacent particles to stick 
together to form large agglomerates on the grain burning 
surface, prior to passing into the flow. At higher burning 
rates (pressures) the particles would be ejected from the 


surface more rapidly, resulting in smaller agglomerates. 
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3. Near=-Propellant Surface Data 

The test results obtained near the grain surface are 
summarized in Table 4.3. D35 values are plotted in Figure 
4.4. As can be seen directly from these results, there was 
a strong dependence on chamber pressure of the meeistietats 
Size produced from the Zrc. At the lower pressures (159 and 
188 psig), the distribution of particulates ejected from the 
grain surface were trimodal and bimodal, respectively. 
Figures 4.5 and 4.6 detail the characterization of the 
particulate distributions. The dominant modes were centered 
at approximately 30-37 microns and 65 microns. 

The Sauter mean diameters determined from the 
scattering data were 63.5 microns for the run at 159 psig 
and 72.8 microns for the run at 188 psig. This was 
considered to be good repeatability, run-to-run, allowing 
fOmmevardatvons in the Combustion process from grain to 
grain. As discussed above, many of the ZrC particles, upon 
reaching the burning surface, apparently react to form Zr0>5 
and agglomerate before being ejected into the gas flow. 

Test number three of the near grain results (248 
psig) showed a reduction in Sauter mean diameter, with a 
measured D3 of 32.2 microns. The Jdestribution of 
particulates at this pressure was still essentially bimodal 
as depicted by Figure 4.7. However, the dominant modes of 
the distribution were now centered at approximately 12-13 


microns and 27 microns. There was also a small number of 
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Particulates in the range of 3-4 microns, indicating that 
some of the particulates at this pressure (248 psig) were 
ejected directly into the gas stream. The 448 psig test 
resulted in a D3, of 15.2 microns. Figure 4.8 reveals that 
the distribution of particulates continued to be bimodal. 
It is apparent from Figure 4.8 that a much larger portion of 
the original metal (ZrC) particles was ejected directly into 
the gas flow, as evidenced by the well-defined mode centered 
between 3-4 microns. The general observation here was that 
as motor chamber pressure waS increased, the size of 
agglomerates on the grain surface decreased, similar to the 
behavior of aluminized propellants reported in References 21 
ama 22. It was also observed that a greater number of the 
Original ZrC particulates were apparently immediately 
ejected into the gas flow as they reached the burning 
surface as chamber pressure was increased. In all cases, 
the distributions of particulates near the burning surface 
were bimodal or trimodal. The quantities of particles in 
the smallest modes of the trimodal data were quite small. 
The most difficult obstacles to overcome in 
obtaining near grain light scattering measurements were: 
(1) keeping the window surfaces clean, as the ignition 
process (due to its close proximity to the windows) made it 
@ifficult to effectively purge these surfaces, and (2) 
multiple scattering effects due to dense concentrations of 


particulate matter in the near grain region. The problem of 
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multiple scattering has been addressed above. The empirical 


corrections presented by Gulder [Ref. 28 ] are for a Malvern 
system operated in the model independent mode. This mode was 
also used in the current investigation. When obscuration is 


approximately 0.99, the measured D2. 1S approximately a factor 
of two too small. Gomi [Ref. 27 } has also discussed the effects 
of multiple scattering. For monodispersed spheres Netzer and 
Powers [Ref. 1 ] have also presented data for correcting D..5 


mem Wagqh Obscurations. 

The more difficult problem that had to be overcome 
in the near grain region was effectively purging the window 
surfaces. The only alternative here was to repeat the 
experiment until usable data were obtained. As noted 
earlier, three of the tests conducted in the near grain 
region were invalidated due to the fact that the window 
surfaces had become contaminated as a result of the ignition 
process. 

4. Nozzle Entrance Data 

Table 4.4 presents the results from the measurements 
obtained prior to the nozzle entrance. These data allow one 
to examine the change in particle size across the nozzle 
when compared to the exhaust data which are contained in 
Tables 4.5 and 4.6. The D359 data obtained using the MALVERN 
are also plotted in Figure 4.16. Figure 4.17 plots these 


same data together with the MALVERN data obtained near the 
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burning surface. Figures 4.9, 4.10 and 4.11 characterize 
the volume distributions associated with the 150 psig class 
of experiments. Note in Figures 4.9 and 4.10 that the dis- 
tributions remained bimodal in nature, with the dominant 
modes being centered at approximately 35 and 68-70 microns 
in the case of Figure 4.9 and at approximately 27 and again 
at 68-70 microns in the case of Figure 4.10. In the case 
of Figure 4.11, note that the profile was trimodal, but 
shifted toward the smaller end, with dominant modes centered 
at approximately 3-4 microns, 20 microns and at 38-40 mi- 
crons. The Sauter mean diameter obtained from the scatter- 
ing data was 30.7 microns (Test No. 3). A post-test 
alignment check done with the MALVERN system in the align 


mode revealed that some significant window contamination had 


occurred, but physical examination of the windows showed 


very little evidence of contamination. The results of this 
one test, although questionable, are included £Otr 
completeness. Comparison of the Sauter mean diameters and 


volume distribution profiles for near grain and nozzle 
entrance measurements completed for the 150 psig class 
reveal only slight differences. This occurred even though 
the residence time (tye.) was approximately 45 milliseconds. 

As discussed above the equilibrium adiabatic flame 
temperature was 2916°K. The actual temperature may have 
been somewhat less than this due primarily to heat loss 


through the uninsulated, stainless steel walls of the 
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combustor, and to the nitrogen gas used to purge the viewing 
Windows. The observed particle behavior would be expected 
if the ZrC is completely consumed near (or on) the 
propellant surface, resulting in liquid Zr0O5 particles 
flowing in a low velocity (M << .1) environment. | 

The NPS scattering data included in Table 4.4 will 
be discussed later in conjunction with the comparison of the 
MALVERN and NPS measurement systems. 

Figures 4.12 and 4.13 represent the volume distribu- 
tions associated with the runs completed for the 250 psig 
class. This represents data for tests 8 and 9 of Table 4.4. 
The Sauter mean diameter (SMD) determined from scattering 
measurements was 49.0 microns corresponding to test number 8 
and 48.0 microns associated with test number 9. This again 
represented good test-to-test repeatability. The dominant 
modes were centered at approximately the same locations for 
both tests. Comparison of these data with that obtained in 
the near grain surface region showed that there was approxi- 
mately a 16 micron particle growth over the length of the 
motor. Note that only one measurement was taken in the near 
grain region. If the point was valid, then particle 
collisions and/or wall collisions/shedding may be suspected. 

Test numbers 10 and 11 represent data obtained for 
the 400 psig class runs. The SMD's determined were 32.0 and 
27.6 microns, respectively. Figures 4.14 and 4.15 are the 


volume distributien profiles for tests 10 and lle 
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respectively. Note that each was monomoda Wein ieee sme 
comparing these profiles to those obtained near the grain 
surface, it was observed that the 3-4 micron diameter 
particles were eliminated as they passed along the length 
of the Combustor. Again, if ZFC 1s assumed to Yveaer to ZrO. 
on or near the propellant surface, particle collisions could 
explain this observation. There was an appreciable increase 
in particle diameter from near the grain surface to the 
nozzle entrance section. Figure 4.16 is a graphical depiction 
of SMD plotted against chamber pressure. Note that the SMD 
decreases with increasing chamber pressure. Figure 4.17 
presents the data for both near-grain and nozzle-entrance 
measurements. From this figure, it is clear that particle 
growth/agglomeration occurred across the length of the 
combustor for pressures above 200 psi. 
>. SEemaust Nozzle Data 

Table 4.5 is a listing of test results obtained in 
the exhaust for the condition of nozzle exit plane pressure 
equal to the critical pressure at the nozzle throat. The 
nozzles used in these tests had a 45° converging section and 
a short flat (.01"-.015") in the throat section. All nozzles 
were truncated at the throat after the .01"-.015" flat. 
Figures 4.18, 4.19, 4.20, 4.21 and 4222 presemesrnesvolum— 
distributions associated with tests one through five, respect— 


Viele Figure 4.23 presents the results obtained on 
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test 6 using the NPS system. The distance downstream of the 
nozzle exit plane varied slightly as noted in Table 4.5 due 
to the geometry of the apparatus and set-up. All the 
distributions were monomodal with the exception of the test 
conducted at the lowest pressure (147 psig). Here a second 
broadly distributed mode was observed, which covered a range 


of approximately 4.5-12.0 microns. Except for this peime 





the SMD varied only slightly as chamber pressure was increased. 


The Weber number is defined as 


We = dag Og (Ug = Wag) 2/9 


where: 


dag = the diameter of an agglomerate (m) 
Pg = the density of the gas 

Uag = the velocity of the agglomerate 
Ug = the velocity of the gas 


0 = the sSumEiace mecns on. 


It is the ratio of inertial forces to surface tension 
forces. Gany [Ref. 22] found the Weber number to be a good 
correlation parameter in the analysis of the break-up of 
aluminum oxide agglomerates in the nozzle environment. Gany 
[Ref. 22] notes that for low Weber numbers the droplets are 
spherical, and tend to distort when the Weber number exceeds 


about four (for aluminum oxide). The distortion continues 
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to increase as Weber number increases until break-up occurs 
for Weber numbers in the range of 12-20 [Ref. 22}. 

Figure 4.24 is a plot of SMD versus motor chamber 
pressure. It 1s apparent from the curve that the 
agglomerates broke up into nearly the same size in the 
converging section of the nozzle, irrespective of chamber 
pressure (unless very low) or particle size entering the nozzle. 

The equilibrium adiabatic combustion temperature was 
approximately the same as the melting temperature for Zr0O>. 
Thus, even though some thermal lag can be expected, the 
particles could be expected to rapidly solidify as the gas 
expansion in the nozzle occurs. This would prevent any 
further reduction in particle size. 

Table 4.6 contains the results obtained in the 
exhaust for the condition of nozzle exit pressure equal to 
ambient pressure. The nozzles all had converging half 
angles of 45° and diverging half angles of 159. Nozzle 
throat diameter and exit diameter were determined using 
lsentropic relations for a ratio of specific heats equal to 
1.2. Figures 4.25 to 4.32 present the volume distributions 
associated with tests 1 through 8, respectively. It is 
interesting to note that all the distributions were 
multimodal with the exception of test number 5, which was 
PUnme-at 215 psig. The primary factor which differentiated 
test five from the other tests was the distance behind the 


nozzle exit plane at which the measurements were obtained. 
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The measurements were all made 2.7 inches downstream of the 
nozzle exhaust plane except for test five, where the 
distance was 0.6 inches. The tests run with the 2.7 inch 
distance were necessary because laser diffraction data were 
being obtained simultaneously through the motor, and the 
physical dimensions of the apparatus limited the minimum 
distance to 2.7 inches. Test 5 results indicated that the 
majority of the particles remained monomodal and essentially 
unchanged in mean diameter (19 microns) from the throat to 
the exit plane. The measurements made further aft had one 
mode centered at the same 19 micron diameter. However, the 
number of these particles were less, and significant numbers 
of smaller particles were observed in a bi- or tri-modal 
distribution. The reduced particle size could be due to particle 
breakup. However, it may have resulted from the reduced 
number of particles in the scattering volume, which was locat- 
ed through the center of the spreading exhaust jet. 

The presence of strong thermal gradients along the 
boundary of the exhaust jet initially caused some problem in 
obtaining exhaust measurements. During “ini¢gial “runes 
scattering data indicated the presence of abnormally large 
particulates in the exhaust, on the order of 127 microns. 
Knowing that particle size entering the nozzle was 
substantially less than 127 microns, typically around /70 
microns for 150 psig runs, a search for causal factors was 


pursued. Conversation with Hughes and Taylor at MALVERN 
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Instruments confirmed the idea that thermal gradients or 
density gradients were indeed biasing the scattered light 
data, i.e., the presence of strong thermal gradients along 
the boundary of the hot gas jet, or density gradients 
associated with a shock. The solution to this difficulty 
was a rather simple procedure. When unstable refractive 
index of the suspending medium is present as a result of 
strong thermal gradients, a steering effect occurs that 
steers light to low angles, hence biasing particle size to 
the large end of the scale. The procedure employed when 
using the Master Sizer program was to inhibit data from the 
affected inner rings, using the "KILLDATA" command. Once 
the user had examined the raw data and observed which rings 
had been affected, then those rings were inhibited by the 
software command "KILLDATA." Care must be exercised in this 
procedure to ensure that no scattering data will be lost as 
a result of inhibiting data from the affected inner rings. 
In the latter case there is no recourse, scattering data at 
that point becomes no longer separable from the steered 
mignt. 

Figure 4.33 is a graphical depiction of particle 
size (SMD) versus motor chamber pressure. There waS a 
SieMewreance1on In pamtiele size with increasing pressure. 
Figure 4.34 plots both exhaust conditions. The reduced size 


in the downstream, spreading plume is apparent. 
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Exhaust measurements were made for both Pe= ce and P,= 


P , to determine if any change in particle size was 


ehainealker\ 
incurred as a result Of the =supersonilcyexpansiony co. seneumiei 
Particle breakup downstream of the nozzle exhaust was observed. 
The limited data at approximately 200 psi indicated that little 
breakup occurred within the nozzle, downstream of the throat. 
In aluminized propellants, one might reasonably expect that 
vapor condensation occurring during the expansion process might 
well affect the size distribution of exhaust particles. How- 
ever, in the case of Zirconium Carbide, the temperature ineere 
chamber is not sufficiently high for appreciable amounts of 
either Zn, 7e@iern ZrO. to be present in a vapor phase. 
6. SEM Analysis 

As a means of independently verifying the scattering 
results, residue was collected from the face of the nozzle 
converging section. This residue was then placed in a solution 
of lab grade acetone and put into an ultrasonic cleaner, mixed, 
and allowed to settle. After repeated cleanings, the acetone 
was evaporated off in an oven under reduced pressure. Using 
this powder, seven samples were prepared for examination with 
the scanning electron microscope. Three of the samples were 
prepared using a non-conducting epoxy base which had to have 
a gold flashing applied in order to provide for a conducting 
surface. The remaining samples were mounted using a conduct- 
ing colloidal carbon. Figures 4.35 and 4. 36nare S35!) pnoces 
graphs taken of the samples. Particles were counted from each 


sample in an attempt to get a representative size distribution. Same 218 


80 


Tn Ue CCL 


PERSONAE ems 
7 i Rw + a — - me . = @ 3 ee 
was rs i ae s ma? = a pe 7+ ot «.. aie qe: Ge, ~~ 


at oo a= 
= = 





Figure 4.35 SEM Photograph of Residue Collected 
from Exhaust Nozzle 
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Figure 4.36 SEM Photograph of Residue Collected 
from Exhaust Nozzle 
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particles were counted with sizes ranging from 1-20 microns. 
The smaller particles (3 microns and less) were difficult to 
see, but did exist in the sample. No particle larger than 
20 microns was observed, but this by no means implies that 
no larger particles were present. Collection of residues 
and/or exhaust products 1S a crude method. When collected 
in small quantities from wall surfaces it cannot yield a 
valid statistical sample. However, it does provide some 
reassuring insight as to what was physically present in the 
sample. 

The particles were grouped strongly in 2 bands, one 
band from 5-12 microns and the other from 14-17 microns. 
The sample was taken from the 427 psig run where the mean 
particle size (SMD) prior to the nozzle entrance was 27.6 
microns and in the nozzle throat region it was 19 
microns. Thus, the SEM data were in reasonably good 
agreement with the optical data. In Figures 4.35 and 4.36 
the larger irregular-shaped particles have reasonable 
sources, one being residue from the binder/inhibitor, and 
the second being residue from the lab grade acetone used to 
clean the samples. 

D. COMPARATIVE EVALUATION OF THE MALVERN 2600 AND THE 
LOCALLY DESIGNED LASER DIFFRACTION SYSTEM 
From the outset of this investigation, a comparative 
analysis of the two diffraction instruments was a primary 


goal. Some of the limitations of the NPS systsm were 
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obvious, based on a review of the system hardware and 
software. 

The differences in software and in data reduction 
algorithms were significant. The NPS algorithm is based on 
a ratioing technique developed by Buchele [Ref. 10). This 
method yields the size parameter D3 by the process of 
fitting a theoretical curve for monomodal distributions to a 
measured intensity profile. The technique is discussed in 
detail in Chapters II and III. The MALVERN Master-Sizer 
operating system, by comparison, inverts the scattered data 
con Obtaln a size distribution. The NPS system could be 
upgraded in this area by adopting one of the quadrature 
methods described by Koo [Ref. 16]. 

The differences in hardware are highlighted by the 
design of the detector system photodiode arrays. The 
MALVERN systems annular array design allows ffor’ the 
detection of much weaker intensities. The field of view of 
the NPS system was also limited by the array dimensions. 

A comparison was made using the NPS laser system to 
measure particle size at the nozzle entrance. Figure 4.37 
1s a plot of comparable through-motor measurements made by 
the MALVERN system and the NPS system in the motor section 
just prior to the nozzle entrance. AS can be seen from the 
plot, the NPS system underestimated the mean particle sizes 
at low pressures where the distributions were multimodal. 


As discussed above, the NPS data reduction algorithm assumes 
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a monomodal distribution. The larger particles in the bimodal 
distribution (60-70 microns), scattered light at angles below 
the range of the NPS detector. To further verify this effect, 
the MALVERN data for the run at 161 psig was reduced, after 
inhibiting data that were outside the field of view of the NPS 
system. Figure 4.38 depicts the resulting distribution. What 
occurred was that an artificial mode appeared at 24.7 microns. 
Figure 4.39 1S a measurement made at a similar pressure (175 
psig) with the NPS system. A SMD of 24 microns was obtained. 
Thus, when distributions are multimodal, accurate values of 
D5 are not obtained with present data reduction algorithm 
used with the NPS apparatus. Figures 4.40, 4.41 and 4.42 are 


included for completeness. These figures represent the data 
used to construct Figure 4.37. At higher chamber pressures 
(400 psig class), where the distributions became more mono- 


modal and the dominant scattered light was within the field 
of view of the NPS system, better agreement of data existed 
(Figure 4.37). 

The conclusion is that the NPS system functions well when 
the distribution of particulates is monomodal and the linear 
Gdiede amecay 1s placed such that the majority of the scattered 
hohe tuwkehametme £ield of view of thesdetector. 

E. SOURCES OF MEASUREMENT ERROR (UNCERTAINTY) 
In the course of this investigation a number of 


measurement errors were encountered from a variety of 
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sources. Measurements were made in both the motor chamber 
and exhaust environments, and each environment presented a 
slightly different set of difficulties. 
Measurements made in the motor environment were subject 
to the following sources of uncertainty: 
(a) scattering by the aperture (window frame) 
(b) limited field of view 


(c) multiple scattering effects due Be dense 
concentrations 


(d) surface reflections from the windows in the motor 


(e) diameter of aperture limited scattering angle to a 
maximum of 6° 


(f) cleanliness of window surfaces 

(g) extraneous light 

(h) quenching effects due to chamber walls, and the 
introduction of nitrogen used as a purge for the 
window surfaces 

(He Inetcting "or Scattered light, 1.@€°7" sample not 
completely within the range of the focal length of 
the receiving lens. In the event vignetting occurs, 


particle size is biased to the larger sizes. 


Scattering by the aperture is discussed in some detail 


by Buchele [Ref. 10]. Although some scattering from this 
AeuLce does Occur, if the particle concentration is 
sufficiently high (as in the present case), particle 


scattering will dominate. 

The limited field of view of the NPS system did 
Significantly impact the scattering data obtained. In the 
case of the present investigation, the physical dimensions 


of the photodiode array did not allow for scattered light to 


ou 


be detected at angles less than .01 radians. Also in the 
current configuration of the NPS scattering device, the largest 
particle size which could reasonably be detected was 40 microns. 
The problem may be overcome in large part by increasing the 
dimension of the array. A 4096 diode RETICON array has been 
purchased, and should be installed in the near future. This 
will provide a much enhanced field of view for the NPS instru- 
ment. 

Multiple scattering effects tend to occur when the separat- 
ion distance between particles becomes small. Correction 
methods have appeared in the literature [Refs. 27,28]. These 
effects were present in the motor measurements and empirical 
corrections of the data were required. A reduced combustion 
chamber width is required to reduce this problem. 

The problems associated with windowed surfaces are three- 
fold in nature. The first is that some reflections are present. 
The MALVERN system in the align mode allows sufficient adjust- 
ment of the detector in order to eliminate the effects due to 
reflections. The second problem was a result of the physical 
diameter of the window which could pass light. This restricted 
the maximum scattering angle to 6°. If very small particles 
had been present, this would have had a biasing effect on the 
received scattering data. The apparent size would have been 
larger than the actual size. The third problem associated with 


the windowed surface was to provide sufficient purge to keep these 


OZ 


surfaces clean. A trial and error method of adjusting purge 
flow rates was employed. If the windows had become 
contaminated in a previous run, the flow rate was increased 
until the surfaces remained clean and usable data were 
obtained. 

The introduction of the purge nitrogen, along with the 
walls of the motor chamber, had some undetermined quenching 
effect on the gas and particle temperatures. Although the 
percent of gas flow comprised by the nitrogen purge was 
estimated to be between 3 and 9 percent for the given range 
of pressures, it was not possible to determine what 
influence the quenching effect had on particle size, or the 
extent to which it actually occurred. 

Extraneous light can also be a factor. Buchele [Ref. 
10] details these effects. The ratioing method adopted in 
the data reduction algorithm for the NPS system has the 
desirable feature of eliminating this effect. Reference 10 
points out the fact that extraneous light is most likely to 
be a problem at small angles. 

Vignetting effects occur when scattered light is not 
captured by the receiving lens. This tends to bias particle 
size toward the larger size bands. The way to overcome this 
difficulty is to ensure that the field of particulates is 
well within the focal length of the receiving lens, so that 


all of the light scattered by the sample will be captured by 


03 


the (Fourier transform) lens and focused onto the detector 
abraye 

The measurement errors that were found to be associated 
with obtaining scattering data in the exhaust are listed 
below: 


(a) unstable refractive index of the suspending medium 
which result from thermal or density gradients 


(b) turbulence along the jet boundary 


(c) shocks associated with a jet which is not ideally 
expanded 


(ad) the location at which the measurements are taken 
ayiate of the nozzle can affect the measured 
aLS tr Powe Lom 


(e) vignetting, similar to that described previously 


(f) extraneous’ light (also discussed previously in 
conjunction with motor measurement errors). 


Thermal gradients play a significant role in measuring 
particle size in hot gas jets. Close to the exit plane of 
the exhaust, the jet was slightly smaller than the MALVERN 
laser beam diameter. As discussed above, scattering data 
were influenced by the thermal gradients along the jet 
boundary. Also as discussed earlier, the "KILLDATA" command 
can be used to inhibit data from the inner rings where the 
steering effect is most prominent. Buchele [Ref. 10] points 
out that turbulence compounds these effects. It was 
impossible to separate the effects of turbulence and thermal 
gradients on the beam steering. In either case, so long as 
the user is reasonably assured that no scattering data are 


present as a result of large particles in the exhaust flow, 
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the "KILLDATA" command may be used to inhibit data from the 
light steered to low angles. 

The presence of shocks in the exhaust jet which is not 
ideally expanded, is to be expected. The density gradients 
in the region of a shock can also adversely affect scattered 
data. If shocks are identified as a problem, a schlieren 
may be used to establish shock location so that on 
subsequent runs the apparatus can be positioned upstream of 
the shock. 

The location at which diffraction measurements are made 
in a supersonic exhaust can influence the resulting 
distribution of particulates which are observed. The 
further back in the jet, the less concentrated the particle 
field becomes. This may have biased the measurements to 
the smaller particles which were present. 

The effects of extraneous light and vignetting have been 
previously discussed and will not be discussed further, 
except to say that the measurement and subtracting of 
background light plays a role in reducing the effect of 


extraneous light. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


The results of this investigation have demonstrated that 
laser diffraction methods can be an effective means for 
determining particle sizes in  solid-propellant rocket 
motors. Measurements were successfully made in both the 
motor chamber and exhaust environments. Each environment 
presented a different set of difficulties. 

The MALVERN and the locally designed diffraction system 
produced data which correlated quite well when the particle 
distributions were monomodal, and the dominant scattered 
light was within the field of view of the NPS linear diode 
array. When the distributions were multimodal, the NPS 
system did not accurately measure the Sauter mean diameter. 

This problem was the result of both a limited field of 
view and the data reduction algorithm. The latter was based 
upon monomodal distributions. The system could be 
Significantly improved by using larger arrays and by using 
data reduction techniques which are based on quadrature 
methods [Ref. 16]. 

The NPS system will always be limited to higher number 
density flows since the linear array only intercepts a very 
small fraction of the light scattered at a given angle. 

The MALVERN 2600 has proven to be well suited for use in 


solid propellant combustion experiments. The beam diameter, 
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Ne Ves larger theaneoprimum fommene Current application. 
The acquisition of a new beam expander to reduce the diameter 
of the collimated beam would enhance both in-situ motor meas- 
urements and exhaust meaSurements. 
Preijo l= lOnurroleciewaxisyimetric motor £€O a two-dimensional 


motor would improve through-motor measurements in the following 


ways: 
(a) advantage of reduced optical depth 
(eeseauceton Or elimination of multiple scattering effects 
that may be present in the motor chamber environment 
due to dense particle concentrations 
(c) increased range of scattering angles 
(d) allows operation to higher pressures. 
The motor data suggested that ZrC reacts to form 2x0, on 
Or near the propellant surface. Surface agglomerations occur- 


red and were obServed to decrease in size as chamber pressure 
increased. 

At the low Mach numbers in the present motor the ZrO. 
particles appeared to remain practically unchanged as they 
passed down the length of the motor port at very low pressures. 
At higher pressures (above 250 psig), particle growth/agglom- 
eration occurred. 

At low pressures the ZrO. particles were multimodal. At 
higher pressures ( greater than 350 psi) the distributions 


approached monomodal. 


97 


ZrO. particles broke up in the nozzle converging section 


Zz 
to a monomodal distribution with a D5 between 15 and 19 microns; 
independent of the particle size entering the nozzle. Limited 
data indicated that no appreciable change in particle size 
occurred in the supersonic nozzle flow. Are ot themmezz ley 
in the expanding plume, the particle size was significantly 
Smaller (2-8 microns). This could be due to particle breakup, 
or to measurement location biasing. 

SEM analysis of particles collected form the converging 


nozzle wall showed particles with sizes in good agreement 


with the optically obtained data. 
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